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Abstract. We investigate the relation between the Doppler factor determined 
from variations in total flux at 22 and 37 GHz, and the apparent transverse ve- 
locity determined from VLBA observations at 2 cm. The data are consistent 
with the relativistic beaming theory for compact radio sources, in that the dis- 
tribution of /Sapp/t^var, for 30 quasars, is roughly consistent with a Monte Carlo 
simulation. The intrinsic temperature appears to be ~ 2 x 10^" K, close to the 
"equipartition value" calculated by Readhead (1994). We deduce the distribu- 
tion of Lorentz factors for a group of 48 sources; the values range up to about 
7 = 40. 



1. Introduction 

Relativistic effects are commonly invoked to explain the superluminal motion 
and high brightness temperature (Tb) seen in compact radio sources. Tb can 
be found from VLBI measurements or from variations in flux density; but these 
two methods depend differently on the Doppler factor. By measuring Tb in 
both ways, the Doppler factor and the Lorentz factor can be deduced, as well 
as Tb(int), the intrinsic temperature in the synchrotron source. Lahteenmaki 
et al. (1999a) have done this for several sets of sources, and found that most 
of the Tb(int) have a range around 10^^ K. This is close to the "equipartition" 
value ^ 5 X 10^° K suggested by Readhead (1994), and the diamagnetic limit 
~ 3 X 10^^ K calculated by Singal (1986). It is at the bottom end of the range sug- 
gested by Kellermann &: Pauliny-Toth for the limit due to the inverse-Compton 
catastrophe, 1 — 10 x 10^^ K. (See also Kellermann, these proceedings, page 185.) 

In this paper we compare superluminal velocities, /Japp, with Jvar, the Doppler 
factor derived from variability. If a simple relativistic beaming theory is correct. 
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there will be a close connection between these two quantities. The new large 
and reliable sample of superluminal sources found in the VLBA 2 cm survey 
(Zensus et al., these proceedings, page 27) allows this to be done with some 
confidence, but see Section 6 for a discussion of the Jvar- We work with veloc- 
ities from the VLBA 2 cm survey and Doppler factors from the Metsahovi flux 
density monitoring program at 22 and 37 GHz (Lahteenmaki & Valtaoja 1999b, 
hereafter LV99). We convert the Metsahovi values to our assumed cosmology, 
Ho = 65kms^^ Mpc~^, ilm = 0.3, Ha = 0.7, but otherwise use the values found 
in their paper. 



2. The Velocities 

We only use velocities from the 2 cm survey which have three or more epochs 
of observation with an appropriately small error to the fit, and which satisfy 
several criteria for morphology. (Kellermann et al. 2003, in preparation.) The 
velocities are determined by a linear regression through the locations of the 
components being followed, and the errors shown in Figure 1 are the curve- 
fitting errors. We are using the fastest component for each source, on the grounds 
that their velocities should be representative of the true flow velocities. Slower- 
moving components, especially those at a bend in the jet, may be dominated 
by backward shock waves. Forward shock waves might also exist, and trying to 
understand their role is also a goal of the survey. Other geometries have been 
suggested for the jet, including a fast "spine" which we would preferentially 
see, surrounded by a slower shell. In this case, we think, the spine would also 
control the flux variations, so that using the fastest (spine) velocity for /3app 
is appropriate. In the next sections we assume that the pattern velocities are 
identical to the flow velocities; we comment on this assumption in Section 7. 

Five sources in common to the (5var and /3app lists were excluded on the 
grounds that they only had components situated at a bend in the jet. We 
believe that in these cases we see a standing shock wave, or perhaps a stationary 
location in a helical jet where the flow is closest to the LOS, and hence boosted 
most strongly. In cither case the measured velocity is a poor indicator of the 
flow velocity, and not useful in looking for relativistic effects. The final sample 
we use contains 48 sources: 4 galaxies, 14 BL Lacs, and 30 quasars. 



3. The Doppler Factors 

The 22 and 37 GHz light curves give the flux density in an outburst and the 
time constant At. As a model we take cAt to be the radius of an optically thick 
sphere. A solid angle is calculated and then the brightness temperature Tb. 
(See Lahteenmaki et al. 1999a for details.) A further assumption is now needed, 
that the measured brightness temperature may be relativistically boosted above 
an intrinsic temperature Tb(int). LV99 choose Tb(int)= 5 x 10^*^ K, from the 
estimates of equilibrium temperature for a self-absorbed synchrotron source cal- 
culated by Readhead (1994). The variability Doppler factor is calculated as 
<Jvar =[rb(var)/Tb(int)]V3. 
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A more direct way to get the brightness temperature of a compact source is 
to measure its diameter and flux density with VLBI. The Doppler factor can then 
be calculated by reference to an assumed Tb(int). In most cases of interest this 
yields a lower limit, because the cores of the sources are only slightly resolved 
with terrestrial interferometer baselines (Kellermann, these proceedings, page 
185). 
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Figure 1 . (a) Monte Carlo simulation of (/3app , S) for a flux-limited sample 
with N=100, a = —1.25. The diagonal line represents the "I/7 cone", where 
the apparent velocity is maximized for a fixed 7. The semi-circle is the locus 
of points for 7 = 30. (b) As in (a) with (5var calculated for Tint = 4 x 10^ K. 
Circles: Quasars; Crosses: BL Lacs; Triangles: Galaxies. N=48 (c) As in (b) 
with Tint = 2 X 10^° K. (d) As in (b) with Ti„t = 1 x lO" K. 7 = 20 is shown 
as the dashed curve. 



An alternative method of finding a Doppler factor involves the inverse 
Compton effect and the ratio of X-ray to radio flux. This method depends 
on knowing the diameter from VLBI, but is only weakly dependent on the X- 
ray flux, and so is close to the pure radio VLBI method described in the previous 
paragraph. A few cases have been studied in detail and give good results (e.g., 
Unwin et al. 1994). 



4. The {Pa.pp,S) Relation 

Lister & Marscher (1997, hereafter LM97) have made Monte Carlo simulations 

of a flux-limited sample of compact sources drawn from a population that has 
a distribution of Lorentz factor N(7) ~ 7". They assume flux boosting of the 
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form S'~ 5^, which is appropriate for the flat-spectrum core. Figure 1(a) shows a 
reahzation for a = —1.25. (See LM97 for a discussion of the exponents.) In this 
case the Lorentz factors are set such that 7 = 30 is the maximum value in the 
sample, so the points all fall under the 7 = 30 curve. Note that a majority of the 
sources lie under the diagonal line /3app = 5, or I/7, as they should, since 
the probability density for selecting a source peaks near 6 = I/27 (Cohen 1989; 
Vermeulen & Cohen 1994, hereafter VC94; LM97). For a fixed 7, the apparent 
velocity is a maximum on this line. 

In Figures 1(b), (c) and (d) we show values of (/3app) <^var) for our sample, 
for three values of Tb(int). We believe that our data are representative of a 
flux-limited sample (Zensus et al., these proceedings, page 27), and LM97 show 
that a = —1.25 gives a reasonable fit to the apparent-velocity distribution of 
the Caltech-Jodrell fiat spectrum survey (Taylor et al 1996). Hence our points 
should lie within a region bounded by a 7 = const curve, as in Figure 1(a). 
In Figure 1(b) 6 points are off scale on the x-axis, and the values of 5var go to 
100. In Figure 1(d) we show two 7 = const curves, and it appears that, as in 
Figure 1(b), the points will not properly fill any 7=const curve. Figure 1(c) 
comes closest to the simulation, and, given the small number of points (30), 
it perhaps provides an adequate match. These comments about Figure 1 are 
confirmed by calculating the fraction / of sources inside the I/7 cone. For the 
Monte Carlo simulation / = 0.80, and for the observations the values for the 
30 quasars are (for increasing temperature) 0.87, 0.77, and 0.50. On this basis 
1 X 10^^ can be excluded, and 2 x 10^° is somewhat better than 4 x 10^. 

5. The distributions of /3app/<^var and 7 

The product ^9 is a useful quantity because the probability density function for 
9 depends only on 76* for 7^ ^> 1 (see VC94, Fig. 7). Hence the measured distri- 
bution of 76 can be compared with a theoretical distribution, or one generated 
by a simulation, to see how closely the observations match the standard theory. 
In the previous section we found /, the fraction of sources with 7^ < 1, and we 
picked rb(int) = 2 x 10^° K because its value of / is close to the expected value. 
However, it is clear that looking at many values of j6 would be more powerful 
than looking at / alone, as this would test the general distribution throughout 
the region below the 7 = const curves. 

We actually use the ratio /3app/5var, which, for 7^ 3> 1, equals 7^. (A general 
result for a relativistic beam is /3app = (3'yS sin 9, where P = v/c.) We form the 
histogram of /Japp/^var, which counts the points between radii from the origin in 
Figure 1. In principle we can mix galaxies, BL Lacs, and quasars even though 
they have different parent populations, if we assume that a, and the slope of the 
luminosity function, are the same for all the populations. These assumptions 
seem somewhat dubious, however, and we only use the quasars. 

Figure 2(a) shows the distribution of /3app/'^var for the 100 sources in the 
simulation. Figures 2(b), (c), and (d) arc for the 30 quasars, for the three 
values of rb(int). The histogram in Figure 2(d) is substantially flatter than that 
in 2(a) and supports our statement that rb(int)= 1 x 10^^ K is too high. In 
Figure 2(b) the distribution is too narrow, with a large spike in the first bin. 
This corresponds to the high values of c^var for this temperature. In Figure 2(c) 



Variability and Velocity of Superluminal Sources 



181 



0.3 
n/N 
0.15 



0.3 
n/N 
0.15 



T — I — I — I — I — I — I — I — I — I — I — I — I — r 



' (a) 

Simulation 

N=100 




(c) 



Tint=2xl0'° K 
N = 30 



_l I I I I I I I L 



T — I — I — I — I — I — I — I — I — I — I — I — I — r 



(b) 



Tini=4xl0^ K H 

N = 30 



(d) 



T.nt=lxlO^^ K 
N=30 




J I I I I I I I I I I L 



0.3 
n/N 
0.15 



0.3 
n/N 
0.15 



1 2 



/?app/5 



Figure 2. Histograms of Papp/S = jO. (a) Results from a Monte Carlo 
simulation, A'^ = 100. (b), (c), and (d) Calculated for 30 quasars, with three 
different values of Tb{int) used to calculate S^a.!- 



the histogram is a better match than the others, but it still has an excess in the 
first bin. We originally discarded 5 quasars because they had slow or stationary 
components located at a bend, and we believed that in those cases the measured 
velocity was not a good measure of the flow velocity. If we had kept those sources, 
they would have been in the first bin and the fit would be worse. In fact we 
discarded only the most egregious cases, and it is likely that there are other 
sources where we have measured a slow pattern rather than the flow velocity. 
This, we believe, explains at least part of the excess in the first bin. 

The number of sources is too small to allow us to pick a particular "best" 
value for rb(int). However, it seems likely that for many quasars the intrinsic 
temperature in the compact synchrotron source is 1 - 3 X 10^° K. For illustrative 
purposes below we use Tb(int)= 2 x 10^'' K, and we apply it to galaxies and BL 
Lacs as well as quasars. 

We now calculate values of 7, using rb(int)= 2 X 10^° K. However, four 
BL Lacs have negative Pa,pp, with values within 1.5a" of zero. We assume that 
/?app + 2cr is an upper limit to the velocity, and calculate the resulting upper limit 
to the Lorentz factor. The values of 7 found in this way are nearly independent 
of how we pick /3app, as can be seen from Figure 1; the 7 = const curves are 
vertical at /3app = 0, where 7 = 5var/2. We could instead adopt a lower limit 
and assume that the negative sources are in the forward jet, but moving towards 
the center. In this case their flux would be de-boosted, whereas in fact they are 
strong sources. Alternately, we could assume that they are inward-moving in 
the back jet. We reject this because there is plentiful evidence that in these 
compact sources we see only forward jets. 
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Figure 3 shows the distribution of Lorentz factor for the 48 sources, with 
upper Hmits denoted by arrows. The quasars have a broad distribution between 
5 and 25, with a few quasars both higher and lower. The galaxies and most of 
the BL Lacs have 7 < 10. However, the two highest 7 have large error bars 
(the two BL Lacs above the curves in Fig. 1) and they are not reliable. The 
quasars and BL Lacs have a wide distribution of 7, but there are insufficient 
data to see if they might match any of the distributions found in the Monte 
Carlo simulations by LM97. Wc note that some of their distributions for 7 
have a minimum in the first bin, as in Figure 3. It is interesting to note that 
Homan (these proceedings, page 35) has estimated 7 for the source 3C 279, by 
analyzing a change in velocity, and assuming that it is due to a change in 6, 
not 7. He finds 7 ^ 15. Our calculation using 5var and /?app gives 7 = 19, 
for rb(int)= 2 x 10^° K. These methods of finding 7 are largely independent. 
Their good agreement adds to our conclusion that for many sources the intrinsic 
temperature is ~ 1 — 3 x lO^'^ K. 
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Figure 3. Histogram of Lorentz factors, calculated from /3app and Jvar based 
on Tint = 2 X 10^° K. Arrows indicate upper limits. TV = 48 



6. Reliability of (^var 

We have several concerns over the reliability of the (5var values. The first involves 
the procedure used to calculate c^var from a light curve. To partially address this 
we made some independent determinations of (5var using the same procedures as 
Lahteenmaki et al. (1999), and Valtaoja, et al. (1999). Wc used the 8 and 15 GHz 
variability curves from the University of Michigan web site, and the 22 and 
37 GHz data from the Metsahovi web site. We found that the higher-frequencies 
gave better results, because there was less blending between outbursts. Wc did 
not try to investigate the effect of non-overlap of the epochs of the Jyar and /3app 
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measurements; but we have seen that for most sources succcsssivc outbursts 
have a similar velocity. Also, Valtaoja et al. (1999) state that the apparent 
brightness temperature of well-defined flares does not vary much in any one 
source. However, the most important factor appears to be the identification 
of individual outbursts, even at the highest frequencies. This identification is 
somewhat subjective, and apparently we were more conservative and our values 
for (5var were consistently about 30% lower than the values published by LV99. 
In a few cases there were larger discrepancies, apparently caused by our missing 
some large flares not yet posted on the Metsahovi web site (Valtaoja, private 
communication). Thus we think that this systematic effect might limit the 
accuracy of (5var to about 30%. 

The second general concern involves the model used for the source. The 
calculation simply assumes that the time constant of the flare is related to a 
radius by r = cAt, and then that the solid angle of the source is = 7r(r/L')^, 
where D is the metric distance to the source. If we were to assume instead that 
the diameter, not the radius, is given by cAt, then the Ti,(int) would have to be 
increased by a factor of 4 to reproduce the simulation in Figure 1(a). Actually, 
cAt should be regarded as an upper limit to r; this might further increase 
Tb(int). Furthermore, to be more realistic, the sphere should be replaced by a 
shock wave. See e.g., Marscher & Gear (1985). 

7. Pattern and Flow Velocities 

In this work we have assumed that the pattern velocity wc measure with the 
VLBA is the same as the flow velocity of the beam. However it is clear that this 
is not always so. For example, we eliminated several cases where the pattern 
is stationary at a bend in the jet, on the grounds that the bright spot is due 
to a standing shock wave, and its zero motion is not representative of the flow 
velocity. 

VC94 made a simple model to study differing pattern and flow velocities; 
namely, that the flow and the pattern have different 7, and their ratio r = 7p/7b 
has a characteristic value. From the early VLBI data they showed that if r = 1 
in all sources, then 7 cannot be constant, and vice versa. From Figure 3 it is 
quite clear that 7 is not constant in the compact sources. That does not mean 
that r is constant; indeed, since we count stationary sources as having r = 0, 
then r has some distribution also. If the distribution of r mostly has values 
below unity, then in Figure 1 /3app must be increased. The distribution of 7 
will shift to higher values, and its shape may change. It is difficult to predict 
whether the best-matching rb(int) will increase or decrease, as that depends on 
the details of the r distribution. 

8. Conclusions 

We have combined variability Dopplcr factors with superluminal motions, for a 
sample of about 50 radio sources. The distribution of Pupp/^vav is roughly similar 
to that found in a Monte Carlo simulation, provided the intrinsic temperature 
in the synchrotron-emitting medium is Tb(int)~ 2 x 10^*^ K. This is near the 
"equipartition temperature", ~ 5 x lO^'^ K suggested by Readhead (1994). It 
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is below the upper limit based on the inverse-Compton effect, 1 — 10 x 10 
(Kellermann &; Pauliny-Toth 1969; Kellermann, these proceedings, page 185), 
and the diamagnetic hmit calculated by Singal (1986), ~ 3 x 10^^ K. 

The galaxies, and most of the BL Lacs, have 7 < 10, when calculated with 
Tb(int) = 2 X lO^'' K. The quasars have a distribution which is flat between 7 = 5 
and 7 = 25, with only a few quasars above and below these limits. 

We are indebted to H. AUer, M. AUer, M. Kadler, and E. Valtaoja for helpful discussions. 
The NRAO and the VLBA are operated by AUI, under cooperative agreement with 
the NSF. This research has made use of data from the University of Michigan Radio 
Astronomy Observatory which is supported by funds from the University of Michigan. 
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